EUROPHYSICS LETTERS 
Europhys. Lett., (), pp. () 



Diabolical points in the magnetic spectrum of Fes molecules 



O 
O 

o 



X 



Anupam Garg 



Department of Physics and Astronomy, Northwestern University 

H 1 Evanston, IL 60208, USA 

C3 ' 



(received ; accepted 



PACS. 75.10Dg - Crystal-field theory and spin Hamiltonians. 
PACS. 75.50Xx - Molecular Magnets. 
PACS. 03.65Bz - Berry's phase. 

m 
m 

Abstract. - The magnetic molecule Fe8 has been predicted and observed to have a rich pattern 
of degeneracies in its spectrum as an external magnetic field is varied. These degeneracies have 
now been recognized to be diabolical points. This paper analyzes the diabolicity and all essential 
properties of this system using elementary perturbation theory. A variety of arguments is gievn 
to suggest that an earlier semiclassical result for a subset of these points may be exactly true 
g ' for arbitrary spin. 

O 

The molecular cluster [Fe802(OH)i 2 (tacn) 6 ] 8+ (or just Fes for short) has a total spin J = 10 
at low temperatures, and is described to a first approximation by the spin Hamiltonian [jj], |^, ^| 

H = k l J 2 x +k 2 Jl-gn B i-il, (1) 

where k\ > k% > 0, and H is an external magnetic field. Thus the axes x, y, and z are hard, 
medium, and easy, respectively. EPR measurements indicate k\ 0.33 K, k 2 ~ 0.22 K. 

In the absence of any applied magnetic field, the spin of the molecule has degenerate classical 
minima along the ±z directions. Application of a field cants the minima away from ±z, but 
the degeneracy is preserved if H is in the x-y plane. This degeneracy is lifted by quantum 
mechanical tunnelling between the low energy orientations. It is of some interest to calculate 
the tunnel splitting A, since tunnelling plays an important role in the low temperature 
dynamics. A few years ago, without knowledge of the relevance of Eq. (|l]) to Feg, it was 
predicted Q that, for H||x, A would oscillate as a function of H, with perfect zeros at certain 
values, and this effect was explained in terms of interference arising from a Berry phase in 
the spin path integral. These oscillations have now been seen by Wernsdorfer and Sessoli [|| 
using a clever technique which enables Landau-Zener-Stuckelberg (LZS) transitions between 
the levels in question. The underlying value of A can be extracted from the observed LZS 
transition rate. 

In addition to the predicted oscillations, however, Wernsdorfer and Sessoli have also ob- 
served oscillations for certain non-zero values of H z as H x is swept. Villain and Fort |6| have 
noted that if H z is chosen properly, these oscillations also represent perfect degeneracy, i.e., 
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A again vanishes exactly at isolated points in the H x -H z plane [or the full three-dimensional 
(H x , H y , H z ) space]. Thus, all the zeros of A are, in fact, "diabolical points" in the magnetic 
field space. (This coinage is due to Berry and Wilkinson [Q, as the shape of the energy surface 
when plotted against two parameters in the Hamiltonian — H x and H z in our case — is a 
double elliptic cone joined at the vertex, which resembles an Italian toy called the diavolo.) 
Formulas for these points have been found by Villain and Fort, and independently by the 
author || (see below). 

Diabolical points are of interest because of their rarity in real-life physical systems. Indeed, 
the von Neumann- Wigner theorem states that as a single parameter in a Hamiltonian is varied, 
an intersection of two levels is infinitely unlikely, and that level repulsion is the rule. It is useful 
to review the argument behind this theorem. Let the energies of levels in question be E\ and 
E2, which we suppose to be far from all other levels. Under an incremental perturbation V, 
the secular matrix is 



V21 E2 + V22 



(2) 



with V21 = V^- The difference between the eigenvalues of this matrix is given by 

[{Eh. -E 2 + V u - V22) 2 + 4|V 12 | 2 ] 1/2 , (3) 

which vanishes only if 

E 1 + V X1 = E2 + V22, V 12 = V 1 * 2 = 0. (4) 

Hence, for a general Hermitean matrix, three conditions must be satisfied for a degeneracy, 
which in general requires at least three tunable parameters. If the matrix is real and symmetic, 
the number of conditions and tunable parameters is reduced to two jl3j . 

An exception to this rule occurs when the Hamiltonian has some symmetry, when levels 
transforming differently under this symmetry can intersect. For the Fes problem, the intersec- 
tions when H||x or H||z can be understood in terms of symmetry [^), but those with both H x 
and H z non-zero cannot. 

The results reported in Ref. J@] are based on a generalization JIo| of the discrete phase 
integral (or WKB) method jn], |l2| , and are asymptotically accurate as J — > 00. Villain and 
Fort use an approximate version of the same method, with the additional condition k\ — fc 2 <C 
k\. These calculations while involving only elementary methods of analysis, still entail the 
development of considerable calculational machinery, and are quite long. Surprisingly, the full 
global structure of the energy spectrum can be obtained by a much simpler method — text-book 
perturbation theory in fc 2 //ci and the field components H y , H z . This is an extension of an 
earlier calculation by Weigert jl4), who analysed the problem for H y — H z — 0. In particular, 
one can rigorously establish the existence of diabolical points, and find formulas for their 
locations via a series of small calculations. It is hoped that the simplicity of this approach will 
make the subject accessible to a wide readership. 

Before proceeding further, it is useful to develop a scheme for labelling the eigenstates of H. 
Suppose first that H = 0, and fc 2 =k\. The states can then be labelled by the eigenvalue m of 
J Zl and the ground states are m — ±J. If fc 2 is now decreased, states with m differing by an 
even integer will mix. If k\ — fc 2 <SC fc 2 , or J is large, the barrier between m — — J and m — +J is 
large (see fig. 1), tunnelling is negligible, and we can find states |m*) which evolve continuously 
from |m), such that are eigenstates of TL to good approximation. This approximation 

will continue to hold if the field H is turned on, as long as |H| <C H c = 2k\J / g^B- 

The first set of diabolical points lies on the line H y = H z =0, because TC is then invariant 
under a 180° rotation about x. Levels with different parity under this operation can intersect 
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Fig. 1. - Schematic energy level diagram when &2 ~ fci. 



as H x is varied. In particular, the pseudo-ground states m* = ±J are exactly degenerate at a 
sequence of H x values as found in Ref. [[| . Since the symmetry is destroyed if either H y ^= 
or H z ^ 0, so are the intersections, and the points are indeed diabolical. The same is true of 
intersections of levels with m* — ±( J — £), where £ is an integer. 

A similar argument applies when H||z, so another set of diabolical points is expected when 
H x = H y = 0. In terms of fig. 1, the states which are degenerate are no longer symmetrically 
located, and it is possible for say, m* = — J to be degenerate with m* = J — 1. The new 
discovery by Wernsdorfer and Sessoli is that the tunnel splitting between these states also 
oscillates as H x is now varied. As mentioned above, these oscillations are also perfect, and 
the corresponding diabolical points are not associated with any obvious symmetry of TL. (A 
similar situation holds in the spectrum of a particle confined to a two dimensional triangular 
region |Q. Apart from an overall size, which only affects the overall energy scale in a trivial 
manner, a triangle is parametrized by two angles. Two sets of diabolical points arise when 
the triangles are isoceles, but the rest appear when the triangles are scalene with no special 
symmetry.) 

We can thus classify the diabolical points by the m* numbers of the levels which are 
degenerate. Let the state with predominantly negative values of m be labelled by m\ , and the 
other state by m\. We define k = m* + J, and k' = — J- In other words, counting from 
0, the fcth level in the left well is degenerate with state number k' in the right well. When 
k, k' <C J, the semiclassical analysis gives the location of the diabolical point as (H y = 
always) 



h x 



\ZT~A 
J 



J -I- i(fc + fc' + l) 



(5) 



h z « |y( fc - fc ')- (6) 

Here, h = H/H c is a reduced field with H c — 2kiJ/gfiB> A = fe/fci, and £ is an integer. 

Another way to label the degeneracies is to number the levels in order of increasing energy, 
starting with 1 for the lowest level, and then simply give the numbers of the two crossing 
levels. Thus if the lowest two levels are degenerate (k — k' = 0), we will say that levels 1 and 
2 cross, while for k = 0, k' = 1, or k = 1, k' = 0, we would say that levels 2 and 3 cross. This 
labelling is not unique, but we will find it convenient. 

With this background, we now turn to our calculations. Following Weigert Q we regard 
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the ki term in Eq. (1) as the perturbation, along with the y and z components of H. It is 
convenient to divide all energies by fci, and write H = Hjk\ = Ho + Hi, where 

Ho = Jl-2Jh x J x , (7) 
Hi = \J% - J(h-J + + h+J-), (8) 

where J± — J y ± iJ z , h± — h y ± ih z . These notations for J± are unconventional, but they 
are now convenient, as we will take the quantization axis to be x, not z. We will label the 
eigenvalue of J x by n. To zeroth order, the energy of state n is given by 

= n 2 - 2Jh x n, (9) 

Levels n and n' are approximately degenerate if Jh x = (n + n')/2. To see if they are exactly 
degenerate when Hi is included, we find the secular matrix V to an appropriate order in 
perturbation theory, and examine the conditions (Q). We do this for a number of different 
cases. 

Case 1 — levels 1 and 2 cross. Let the degenerate levels be no and no + 1, so that 
Jh x rj (no + For brevity, we label the states by A and B, and denote the matrix elements 
(no + 1| J-\no) etc. by a±, 02, 03, etc., as indicated in fig. 2. Note that all dj can be chosen as 
real. To first order in A and h±, 

V AA = X[J(J + 1) - n 2 ]/2, (10) 
Vbb = \[J(J+1)- {n + l) 2 }/2, (11) 
V AB = -Jh+a 2 . (12) 

The conditions for diabolicity are thus 

Jh x = (n + ±){1-1\), h y = h z = 0. (13) 

Writing no = J — £ — 1, this is identical to Eqs. rtq) and (^) with k = k' = 0, once we recognize 
that (1 - A/2) = (1 - A) 1 / 2 + <3(A 2 ). Since — J < no < J — 1, there are 23 such points. 

The conclusion that these points lie on the line h y = h z = is unchanged if we go to higher 
order. The relevant condition is clearly that for off-diagonal elements. Contributions to the 
AB element of the second order secular matrix arise from intermediate states no + 2 and no — 1. 



ANUPAM GARG: DIABOLICAL POINTS IN THE Fc 8 SYSTEM 



5 




Fig. 3. - Energy level diagram for case 2. 



A short calculation gives V^g = \a2J{a 2 + a§)/i_/8. Adding this to Eq. ( |l2| ) and setting the 
sum to zero, we again obtain the conditions h y — h z = 0. 

Case 2 — levels 2 and 3 cross. Let the lowest energy level be no, and let no ± 1 be 
approximately degenerate. This requires Jh x ~ no. Again, we denote the states no ± 1 by A 
and B, and the various matrix elements of J± by a\ to 04 as in fig. 3. To O(A), Vaa and Vbb 
are given by A[J( J + 1) — (ri ± l) 2 ]/2. The order h 2 , hi contributions to the diagonal terms of 
the second order secular matrix are found to both be equal to 4J 2 (h 2 + hl){J(J+l)~ no + l]/3. 
The interesting terms are Vab and Vba- Including first order pieces from XJ 2 , and second 
order pieces from h±, we get 

Vab = QA+ J 2 h 2 + )a 2 a 3 . (14) 
For a diabolical point, therefore, the vector h must have components 

1 



J 



n (l- ±A),0,±\/A . (15) 



With no = J — i — 1, these are exactly the lowest order terms in an expansion in A of Eqs. 
and @ with k = 1, k' = 0. Since — J + 1 < no < J — 1, there are 2 J — 1 such points. 

Case 3 — levels 3 and 4 cross. This case can arise either with k = k' = 1, or with 
fc = 2, fc' = 0, but we shall be able to distinguish between these. Refering to fig. 2 again, 
the degenerate levels are no — 1 (C) and no + 2 (D). Equality of Eq and Ed again requires 
Jh x « (n + 5). To first order in A, Vcc ~ Vdd = — 3A(n + i), so that the diagonal elements 
are equal when 

Jh x = (no + \){l-\\). (16) 

As in case 2, it is the off-diagonal term which is of greater interest. The secular matrix is now 
diagonal in first order, and off-diagonal terms only arise in second and higher orders. Second 
order terms arise from the combination of one h±J^ term and one XJ 2 term, while third order 
terms arise from three h±J^ terms. The net result is 

V C d = -\(h 2 + J 2 + \)h + Ja ia2 a 3 . (17) 

This can vanish in two ways. The first is to have h y = h z = 0, in which case the diabolical 
field is given by Eq. ( |T3| ) again. This case corresponds to k = k! = 1. 
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The second way for Vcd to vanish is for the factor in parentheses in Eq. (|T^) to vanish. 
This happens when 

h y = 0, h z = VA/J. (18) 



In conjunction with Eq. (16), this is seen to be the same as Eqs. (||) and (||) with k = 2, k' = 0, 
and no = J — ^ — 2. 

It is clear that this procedure gets rapidly more tedious if we apply it to cases with larger k 
and k' . It is more useful to consider higher order perturbative corrections for the cases treated 
above. In the argument leading to Eq. (pf), e.g., we have only gone up to O(A). It is obvious 
that inclusion of higher order terms can at best alter the value of h x at the diabolical point 
by terms of order A 2 , h yl and /i 2 , but cannot destroy the existence of a perfect degeneracy. 
The same argument applies in all the other cases, and constitutes a constructive proof of the 
existence of diabolical points. 

It is particularly interesting to investigate the subset of diabolical points on the line H y = 
H z — in greater depth. As noted before, these points correspond to k = k', and the 
degenerate levels have n quantum numbers differing by an odd integer. Thus they can never 
be coupled by the remaining perturbation TL\ = XJy, and the problem is effectively one of 
non-degenerate perturbation theory. Let us consider case 1 first. The second order correction 
to the energy of state A arises from the intermediate states no ± 2, and to that of state B from 
uq — 1 and uq + 3. It suffices to find the energy denominators assuming that Jh x = (no + 5). 
A short calculation gives 



V ( bI ) 3 



.2 ., m 2 ,oJ 4 f -n - 1 

5n + 2 



[J{J+ 1) - K + n + 1)1 + 2n 2 + 



(19) 



Along with Eqs. (P)-(|rij), this means that to 0(A 2 ), the states are degenerate when 

Jh x = {n + \){\-\\-\\ 2 ), (20) 

which is precisely what Eq. (^) also gives. 

In the same way, for the subcase k = k' = 1 of case 3, we obtain 



<cc j = -_La 2 
V§1 ) 40 



Vnn \ \2 



[J(J + 1) - (n 2 + n - l)] 2 - 6n 2 + ( _\ Q 



(21) 



Including lower order terms, the condition for degeneracy is found to be identical to Eq. (p0|). 

The fact that the two pairs of states k = k' = 0, and k = k' = 1 are simultaneously 
degenerate (at least to order A 2 ), is very striking. Calculations to 0(A 2 ) were in fact done 
by Weigert JL4|, but he did not perform them sufficiently explicitly, and reached the opposite 
conclusion, i.e., that the degeneracy conditions would be different. It is clear, however, that 
this equality is a result of the simple form of H, and is violated when higher anisotropics such 
as ( J x ± i J y ) 4 are included. 

The second striking feature about the result ( pp| ) is that there are no terms like A 2 J 2 or 
A 2 nQ etc. on the right hand side, and that it is agrees precisely with the semiclassical answer. 
Since the latter is obtained in a very different limit, namely, J — > 00, it begins to raise the 
suspicion that it might be exact. To test this suspicion, we have carried the calculation for 
case 1 to order A 3 . For this, not only must we find V^l an d Vg^, but we must also keep 
O(X) corrections in the energy denominators in the calculations for V\j± and Vg^, since Jh x 
depends on A at the diabolical point. The resulting calculation is lengthy, but is efficiently 
done using MAPLE. Almost miraculously, all powers of J multiplying A 3 cancel, as do terms 
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A 3 ?1q with j > 2, and the contribution to Ea — Eb is just A 3 (2no + 1)/16. The condition for 
degeneracy thus becomes 

J^ = (no + i)(l-iA-|A 2 -^A 3 ). (22) 

It will not have escaped the reader that the last factor equals (1 — A) 1 / 2 to 0(A 3 )! 

It is useful to consider the structure of the perturbation series to higher order in A. It is 
clear that we cannot get negative powers of J in the formula for Jh x ; instead it generates 
positive powers. Although the low order analysis suggests otherwise, in principle we should 
expect terms such as X N J K (J + 1) K with < K < N — 1 in iVth order. Such terms would 
be reminiscent of an asymptotic series, and would signal a zero radius of convergence. Such a 
situation would be very odd in our problem since the perturbation A J 2 does not appear to be 
singular. Although plausible, this is far from a complete argument that such terms are in fact 
absent, since we have not excluded terms such as A "n*" 1 in iVth order. 

Further evidence that the result (||) is exact comes from looking at low values of J. We 
have done this for J up to 2. For J = 1/2, there is nothing to prove as the only degeneracy 
is at h x — 0, which is also guaranteed by Kramers's theorem. For J = 1, the energies are 
directly found to be E ±1 = 1 + |A qp (4/i 2 J 2 + \\ 2 ) 1/2 , and E = A, so E\ = E Q when 
Jh x = (1 — A) 1 / 2 /2. For J = 3/2, H separates into two 2x2 matrices in the J x basis, which 
we call Mi and M^- Both eigenvalues of Mi coincide with those of Ma at h x — 0. This is again 
Kramers's degeneracy. In addition, one eigenvalue of Mi coincides with one of M2 precisely 
when h x = 2(1 — A) 1 / 2 /3. For J — 2, 7i separates into a 3 x 3 matrix (Mi) and a 2 x 2 matrix 
(M 2 ). The expected degeneracies are at h x = (1 — A)^ 2 /4 and 3(1 — A) 1 / 2 /4. At the second 
value of h x , one eigenvalue of Mi indeed coincides with one of M 2 . At h x = (1 — A) 1 / 2 /4, 
however, two distinct M\ eigenvalues coincide with two M2 eigenvalues. Thus we again see 
the simultaneous degeneracy of two sets of levels (k = k' = 0, and k = k 1 = 1), leading us to 
believe that this feature is also generally true. A rigorous proof of these conjectures remains 
an open problem. 

This work is supported by the National Science Foundation through Grant No. DMR- 
9616749. 
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